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I. INTRODUCTION 



In nuclear and high energy physics [1,2], it is necessary to obtain the exact bound state 
energy spectrum of non-relativistic or relativistic wave equations for various potentials. In 
relativistic quantum mechanics, the Klein-Gordon and Dirac wave equations are most fre- 
quently used to understand the dynamics of a particle in large class of potentials. Therefore, 
much works have been done to solve these relativistic wave equations for various radial and 
angular potentials. Unfortunately, some quantum mechanical equations with these poten- 
tials can be solved analytically only for s-wave equation (i.e., I — states). Some other works 
have been done to study the I ^ bound state solutions of the Klein-Gordon equation with 
pure scalar, pure vector and mixed vector and scalar potentials in various dimensional space 
using different methods. The scalar coupling constant is almost taken to be equal to the 
vector coupling constant (i.e., 5*0 = Vq) [2-14]. For this case, the Klein-Gordon equation 
reduces to a Schrodinger-like equation and thereby the bound state solutions are easily ob- 
tained by using the well-known methods developed in nonrelativistic quantum mechanics 
[2]. On the other hand, it has been shown that when 5*0 > T^o, there exist real bound state 
solutions (i.e., the radial wave function must satisfy the boundary condition that it becomes 
zero when r ^ oo and also finite at r = 0). However, there are very few exact solvable cases 
[15-19]. The bound state solutions for the last case is obtained for the s-wave Klein-Gordon 
equation with the exponential [15] and some other classes of potentials [16]. Chen et al. [19] 
used the exponential function transformation approach along with an approximation for the 
centrifugal potential to find bound-state solutions of the Klein-Gordon equation with the 
vector and scalar Hulthen potentials for different I. 

The D-dimensional Schrodinger and Klein-Gordon wave equations have been solved for 
various types of angular and radial dependent potentials using the Nikiforov-Uvarov (N-U) 
method [20-25]. The exact bound state solutions for the Schrodinger, Klein-Gordon and 
Salpeter wave equations, in one-dimension {ID), with the generalized Woods-Saxon poten- 
tial are also obtained [26-29]. These solutions of the D-dimensional Klein-Gordon equation 
describing a spin-zero particle are obtained for the case Vq = -S'o ring-shaped Kratzer-type 
and pseudoharmonic potentials [21,22]. Furthermore, analj^ic solution of radial Schrodinger 
equation has been given for the / ^ general Hulthen potential within an approximation to 
the centrifugal barrier term [25,29]. Simsek and Egrifes [30] have investigated the reality of 
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exact bound states of the ID Klein-Gordon equation with complex and/or PT-symmetric 
non-Hermitian exponential-type like general Hulthen potential. Berkdemir [31] has also ap- 
phed the method to solve the Klein-Gordon equation of a spin-zero particle for V{r) = S{r) 
Kratzer-type potentials. Qiang et al. [32] have given an approximate analytic solution of the 
I 7^ Klein-Gordon equation for the general Hulthen-type potential within the approxima- 
tion given in [29] . Saad [33] has presented an approximate solution of the I ^ bound states 
of the Klein-Gordon equation in D-dimensions with the general Hulthen-type potentials fol- 
lowing the model used in [29,32]. Chen et al. [34] have employed two semiclassical methods 
to determine the bound state energy spectrum of Klein-Gordon equation when Sq > Vq. 
Dong ct al [35] have also obtained the bound state solutions of the Schrodinger equation 
for an exponential-type potential following Refs. [29,32]. Very recently, we have obtained 
the bound-state solutions for the D-dimensional Schrodinger equation with the Manning- 
Rosen potential [36] which can be reduced to the Hulthen potential and the exponential-type 
hyperbolic potential [37] using a novel approximation to the centrifugal term [29,36,37]. 

The purpose of the present paper is to study the relativistic characteristics of the scalar 
and vector general Hulthen-type potentials written, respectively, as follows [32,33]: 

= -T^^. Sir) = ro = , ^ (1) 

where Vq and So represent the coupling constants of the vector and scalar general Hulthen- 
type potentials, respectively, a is the screening range parameter, tq represents the spatial 
range and q is the deformation parameter [26-28]. This potential has been widely used 
in a number of areas of physics. In atomic physics, Vq = Ze^/rg where Z is the atomic 
number. Equation (1) behaves like a Coulomb potential Vc{r) — —Ze^/r when r <^ tq, 
but decreases exponentially in the asymptotic region when (r ^ tq), so its capacity for 
bound state is smaller than the Coulomb potential. The Klein-Gordon equation has been 
solved analytically in [30] only for the s-wave. This equation cannot be solved analytically 
for / 7^ because of the centrifugal term, Therefore, for / 7^ 0, we must use an 

approximation for the centrifugal term similar to other authors [29,38-43] to find approximate 
analytical solutions of the D-dimensional Klein-Gordon equation with vector and scalar 
general Hulthen-type potentials. However, the resulting analytic solutions are only valid for 
g = 1 in the I ^ case and for any q value in the I — and D = 1, 3. 

This paper is organized as follows: In Section II, we will derive I ^ state solutions 



within the approximation given in [29] for the D-dimensional Klein-Gordon equation de- 
scribing spin-zero particle with the vector and scalar general Hulthen-type potentials. We 
further separate the wave equation into radial and angular parts. Section III is devoted to 
a brief description of the N-U method. In Section IV, we present general solutions to the 
radial and angular equations in D-dimensions. We also give discussions of the solution in 
various dimensions and the constraints for obtaining the real bound state solutions. Further, 
numerical values for the three lower-lying states when pure vector and pure scalar potentials 
are used with coupling constants Vq = 0.25 and Sq — 0.25 and unity mass particle. Finally, 
the relevant conclusions are given in Section V. 

II. THE KLEIN-GORDON EQUATION WITH SCALAR AND VECTOR POTEN- 
TIALS 

The time-independent Klein-Gordon equation (in any arbitrary £)-dimension) with scalar 
S{r) and vector V{r) potentials, r = |r| ,describing spin-zero particle of rest mass, mo, can 
be written as (in the relativistic natural units h — c — 1) [2,21,22] 

{Vi + [E^i - V{r)f - [mo + S{r)f} = 0, V|, = ^, (2) 

j=l 3 

where Eni denotes the relativistic energy and V|) denotes the D-dimensional Laplacian. 
Usually, it is required that 5*0 > Vq and mo > Eni for the existence of real bound states 
[15-17]. Furthermore, the x is a D-dimensional position vector, with the unit vector along 
X is usually denoted by x = x/r, is the hyperspherical Cartesian components xi, 0:2, • • • ,xd 
given as follows (cf. Refs. [20-25] and the references therein): 

Xi — r cos 9i sin ^2 • • • sin 

X'i — r sin 9i sin ^2 • • • sin 9d-i-i 
X3 — r cos 02 sin ^3 • • • sin 9d-i, 

Xj — r cos 9j-i sin 9j ■ ■ ■ sin 9d-i, 3 < j < D — 1, 
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xd-1 = rcos9D~2sm9D-i, 

D 

XD = rcosdD-i, ^x] = r'^,D>2. (3) 

We have Xi = rcosOi, X2 = rsin^^i, 9i = ip, ior D = 2 and Xi = rcos9i, X2 = rsin^i sin6'2, 
X3 = r cos 92, 92 = 9, for D = 3. The volume element of the configuration space is given by 

D D-l 

Yl dxj = r^-^dhdn, dn^\[ (sin 9jy-^d9j, (4) 
j=i 3=1 

where r e [0, 00), ^1 G [0, 27r] and 9j e [0, tt], j G [2, D-l]. 

The wave function ■i/'j^l^.^'^^''' (^) with a given angular momentum Z can be decomposed 

as a product of a radial wave function Ri{r) and the normalized hyper-spherical harmonics 

Y,%,J9^,92,-- -,00-1)^.8 

^'-TZ'li^) = Rim%^J9„ 92,---, 0D-i), (5) 

where 

Yi^\_, {9i, 92,---, Od-i) = H{9,)Hi92, • • • , 9n-2)H{9n-i), (6) 
which is the simultaneous eigenfunction of : 

L',Y,^!l,,J9,,92r--,9n-i) = lAij+J-mtin-^^^^^ 

Z = 0, 1, • • • , Zfe = 0, 1, • • • , Wi, j e[l,D- 1], A; e [2, L> - 2], 

Ll-,Y,%^_^{9,, ^2, • • • , ^D-i) = Z(Z + - 2)r,f.),^_^(^i, ^2, • • • , ^D-i). (7) 



The angular momentum operators are defined as 
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Ll= y Ll, = ( ^^^'''9k4r) + 2<k<D-l, 

" sin'^-i^fea^ik V 99 J sin^^fe' " " 



' _d__ d 



dXb dxa\ 

Employing the method of separation of variables and substituting Eqs. (5)- (7) into Eq. (2), 
we obtain the following Schrodinger-like equation: 

Ri'i^) + ^^^^K^) + 1 ^^^+^-^) + [E^, _ v{r)f - [mo + S{r)f^ Ri{r) = 0, (9) 



where l{l + D — 2)/r^ is known as the centrifugal term. Furthermore, using Ri{r) — 
7--(^-i)/2g,(^^^^ it is straightforward to find out the radial wave equation: 



g"(r) + [Er^i - V{r)f - [mo + S{r)] 



2 {D + 2l-l){D + 2l-3) 

4^2 



g{r) = 0. (10) 



Obviously, Eq. (10) can be analytically solved [17,27,30] only for / = (s-wave) case. 
For I 7^ case, we have to use an approximation for the centrifugal term similar to the 
non-relativistic cases which is valid ior q — 1 value [29,32,33,38-43]: 



We follow the model used by Qiang et al. [32] and rewrite Eq. (10) for the scalar and vector 
general Hulthen-type potentials (1) as. 



9"{r) + 



2 , 2 {moSo + E^iVo) e""'' {Si - V^) e'^- {D + 21- 1){D + 2/ - 3) 



1 — qe~ 



4^2 



(1 - ge-"'-) 

X g{r) = 0. (12) 

On the other hand, the other angular equations, obtained from the separation procedures, 
are [21,22]: 



sin^-^ Oj dOj 



d9j 



sin^ 9i 



1 



sin^-2 Od-i dOo-i 



I -0-2 d 

sm Hd—i 



+ l{l + D-2)-- 



^2 



sin 9d-i 



del 

It is well-known that the solution of Eq. (15) is 

1 



ll + llH{{e,)^0. 



exp(±i/i^i), /i = 0,1,2, 



H{92,--- ,0d-2) = 0, 
(13) 

H{eD-i)^0, (14) 
(15) 

(16) 



Hence, the above Eqs. (12)-(14), have to be solved by using N-U method [20-31,36] which 
is reviewed briefly in the following Section. 
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III. NIKIFOROV-UVAROV METHOD 



The N-U method is breifly outhned here and the details can be found in [20-31,36,37,44]. 
N-U method is proposed to solve the second-order differential equation of the hypergeometric 
type: 

^^^'^ + ^/-^'^ + ^"^"^"^ = °' ^^^^ 
where a{s) and a{s) are polynomials, at most of second-degree, and t(s) is a first-degree 
polynomial. Using a wave function, ipn{s), of the simple form 

i^nis) ^ (f>n{s)y„{s), (18) 

reduces Eq. (17) into an equation of a hypergeometric type 



a{s)y':,{s) + r{s)y'^{s) + Ay„(s) = 0, (19) 

where 

0'(S) 

t{s) = t{s) + 27r(s), r'(s) < 0, (21) 
and A is a parameter defined as 

A = A, = -nr'is) - I?^I^a"(s), n = 0, 1, 2, • • • . (22) 

The polynomial t{s) with the parameter s and prime factors show the differentials at first 
degree be negative. The other part yn{s) of the wave function Eq. (18) is the hypergeometric- 
type function whose polynomial solutions are given by Rodrigues relation 

TD Jn 

where Bn is the normalization constant and the weight function p{s) can be found by [44] 

u;\s) - ^u;{s) = 0, u;{s) = a{s)p{s). (24) 
The function tt{s) and the parameter A are defined as 



a'[s]—T[s} , I / a'[s] — T[s] ' 



7r(g)^ ' \ ±1 o -a{s) + ka{s), (25) 



A = A; + 7r'(s), (26) 

where 7i{s) has to be a polynomial of degree at most one. To obtain k, the expression under 
the square root sign in Eq. (25) can be arranged to be the square of a polynomial of first 
degree [44]. This is possible only if its discriminant is zero. Finally, the energy eigenvalues 
are obtained by comparing Eqs. (22) and (26). 



IV. SOLUTIONS OF THE RADIAL AND ANGLE-DEPENDENT EQUATIONS 



A. The Z)-dimensional angular equations 

At the beginning, we rewrite Eqs. (13) and (14) representing the angular wave equations 
in the following simple forms [21,22] : 



d0^ 



+ {j - l)ctg9^ 



de, 



+ A 



A 



i-1 



sin'' 9^ 



H{9j)^0, je[2,D-2],D>3, (27) 



d9l_. 



d9 



D-l 



l{l + D-2) 



A 



D-2 



sin^ 9d-i 



H{9d-i)^0, (28) 



where Ap = lp{lp+p — I), p = j — 1, j with Ip = 0,1,2, ■■ ■ are the angular quantum numbers 
in D- dimensions. Employing s = cos9j, we transform Eq. (27) to the associated-Legendre 
equation 



fH{s) 
ds^ 



js dH{s) Aj - Aj_i - Ajs"' 



1 — ds 



+ 



H{s)^0, je[2,D-2],D>3. 



(29) 



By comparing Eqs. (29) and (17), the corresponding polynomials are obtained 

r(s) = -js, a(s) = 1 - s^, a(s) = -Ajs^ + Aj - Aj_i. (30) 

Inserting the above expressions into Eq. (25) and taking ct'(s) = —2s, one obtains the 
following function: 



J 



+ Ai-k 



s^ + k- Aj + Aj_i. 



(31) 



Following the method, the polynomial n{s) is found to have the following four possible 
values: 



7r(s) = < 



^ + Aj_ij s for ki = Aj - Aj_i, 

- Aj-i^ s for ki = Aj - Aj_i, 

+ A,_i for k2 = A,- + {i^Y , 

s - Aj^i for k2 = Aj + , 



(32) 



2 

0-2) 
2 



where Ap — lp-\-{p — l)/2, withp — j — j and j e [2, D — 2], D > 3. Imposing the condition 
t'{s) < 0, where r' = T'(Zj_i), we may select the following physical solutions: 

ki = Aj - Aj_i and 7r(s) = -Ij-is, j e [2, D - 2], D > 3, (33) 

which yields 

t(s) = -2(1+A,_i)s. (34) 

Using Eqs. (22) and (26), give 

A„ . = 2nj{l + Aj_i) + {nj - 1) , (35) 
A = Aj - Aj_i - (36) 

and setting A = A„^., we obtain 

Substituting Eqs. (30), (33) and (34) into Eqs. (20)-(21) and (23)-(24), give 

(t>{s) = (1 - S'f-''^ , p{s) = (1 - s2y.-l+0-2)/2 _ ^3g^ 

Further, substituting the weight function p{s) given in Eq. (38) into the Rodrigues relation 
(23) gives 

y„,(s) = K, (1 - sr'-' §-{1- sT^''- , (39) 
where A^^ is the normalization factor. Finally, the angular wave functions are 

I ■ —n ■ 

H^, (ej) = iV„^ (sin ej)' ' p{lj-nj+ij-2)/2,l,-nj+ij-2)/2) ^^^^ q^^^ ^^q^ 

where the quantum numbers rij are defined in Eq. (37) and the normalization factor is 



^ _(2^^+j_Jhi^ 

Likewise, in solving Eq. (28), we introduce a new variable s = cos9]j_i. Thus, we can also 
rearrange it as the universal associatcd-Lcgcndre differential equation 

d'Hjs) {D-l)s dH{s) u{l-s')-An-2 rrf. . .... 

where 

u^l(l + D-2). (43) 



Equation (42) has been recently solved in 2D and 3D by the N-U method in [45,46]. By 
comparing Eqs. (42) and (17), the corresponding polynomials are obtained 



t(s) = -(D - l)s, (j(s) = l-s^, a{s) ^ -vs^ - Kd- 



(44) 



Inserting the above expressions into Eq. (25) and taking cr'(s) = —2s, one obtains the 
following function: 



\ 



D-3 



+ u-k 



+ k - u + Ad-2, 



(45) 



which gives the possible solutions: 

' + m') s for /ci = i/ - Ad-2, 
(^ - m') s for /ci = 1/ - Ad-2, 



(46) 



CD-3) 



- m' for k2 = I' + (-^) ^ , 



where m' — Id-2 H — j^- Imposing the condition t'{s) < 0, where r' = t'{Id-2), one may 



select: 

ki = p - Ad-2 and tt{s) = -Id-2S, 

which yields 

t{s) = -2(l + m')s. 
The following expressions for A are obtained, respectively. 



(47) 
(48) 



A = Xno-i = 2ni5-i(l + m') + nD-i(ni?_i - 1), (49) 

X = l{l + D-2)-lD-2ilD-2 + D-2), (50) 
and from A — Xno-n obtain the angular momentum quantum number: 

I = lD-2 + nD-l- (51) 

where no-i — n and Zi = m in 3D space [46]. Using Eqs. (20)-(21) and (23)-(24), we obtain 

0(s)^(l-s^)'-^/\ p{s)^{l-sT'- (52) 
The Rodrigues relation (23) gives the following wave functions: 



ynn-iis) = 5nD_i (l - S^) 



2\ 



2\no_i+m 



(53) 
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where Bnjy_^ is the normahation factor. Finally, the angular wave functions are 

Hnu-AOo-i) = N^u-. (sin^z._i)'-- Pi^lf )(cos^i,_i), (54) 
where the normalization factor is 



^ (2n.-. + 2m^ + l)n._,! 
"^-^ Y 2r (nij_i + 2m') ' ^ ' 

with m' is defined after Eq. (46). 

B. Bound States of the Z)-dimensional Radial Equation 

The L>-dimensional Klein-Gordon radial energy eigenvalue equation with the vector and 
scalar general Hulthen-type potentials can be rewritten as 



9"{r) + 



^2 2,2 {moSo + EniVo) e" 



1 - ge-"'' 

\{D + 2l-l){D + 2l- 3)«2e-"^ + {Si - V^) e'^"'' 



^(r), ^(0) = 0, (56) 



(1 - ge-"'^) 

On account of the wave function g{r) satisfying the standard bound state condition, which 
is when r ^ oo, the wave function g{r) 0. Equation (56) can be further transformed by 
using a new variable s — qe~°'^ (r e [0, oo), s e 0]), 

9"is) + ^7^5'(«) + 777^ + (A - 7 + Ki> - iPi + P2 + eiy] 9{s) = 0, 

(57) 

where 

_ y/ml - Ell a _ 2 (^O'S'o + ^mV^) g _ Si - Vj _ {D + 21 - 1){D + 21 - 3) 

^nl — , Pi — 2 ' '^2 — ^ — ) 7 — Z ■ 

(58) 

For bound states, Eni < rriQ, Eni > [17,47,48]. Comparing Eqs. (57) and (17), the 
corresponding polynomials are obtained: 

t{s) = 1 - s, a{s) = s(l - s),a{r) = + (A - 7 + - (A + (^2 + eli)s\ (59) 

The substitution of Eq. (59) into Eq. (25) and taking a'{s) — l — 2s, give the polynomial: 

7r(s) = -| ± 1^[1 + 4(/3i + /32 + el - k)] .2 + 4(/. - A + 7 - 2el,)s + Ael^. (60) 
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If the expression under the square root in Eq. (60) is set equal to zero and solved for k, we 
obtain: 

A; = /3i-7±£„;Vl + 4(/32 + 7)- (61) 
In view of that, we arrive at the following four possible functions of 7r(s) : 



2 ~l~ ^nl 
9 ^nl ~l~ 



2 ~l~ ^nl 



s for A;i - /3i -7 + £„,v^mp7T7), 

s for A;2 = /3i - 7 - £nz a/I + 4(/92 + 7)- 

(62) 

The correct value of 7r(s) is chosen such that the function t(s) in Eq. (21) must have negative 
derivative [44]. So we can select the physical values to be 



£ni - IVTTW^TT) 

iVl + 4(/92 + 7) 

£ni + |Vl + 4(/?2 + 7) 



s for A;i = /5i - 7 + Sni^/l + 4(/52 + 7), 



k = /5i-7-£„iV'l + 4(/52 + 7) and 7r(s) = - 



£n/ + ^Vl + 4(/52 + 7) 



5, (63) 



which yield 



r(s) = 1 + 2£, 



1 + £ni + ^Vl + 4(/92 + 7) 



r s 



l+£ni + ^\/l+4(/92 + 7) 



< 0. 



(64) 



Using Eqs. (22) and (26), the following expressions for A are obtained, respectively. 



X — \n — n + 



1 + 2£„, + Vl + 4(/92 + 7) ^, (n = 0,l,2,---), 



(65) 



A = /9i - 7 - + 2£ni) [1 + Vl + 4(/92 + 7) 



(66) 



where n is the radial quantum number. Solving the last two equations and using (52 — 
7, give 



^nl — 



(/3i - 7 - n^) - {2n + 1)5± _ 4g - - (2n + 1)5±] - (D + 2/ - 1)(D + 2Z - 3) 
2{n + S±) ~ 8q{n + 5±) 

2q{moSo + E^iVo) + Si -V^ n + 5± 



2g2Q,2(n + 5±) 



-, (n = 0,l,2. 



5^^^(l±'^y a^^q^ + ^^^l-^ + q{D + 2l-l){D + 2l-3), (i = 0, 1, 2, • • • ), 

(67) 
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with 5 — 5+ ior q > and 5 — 5- ior q < 0. Furthermore, from setting e^^^ — e^i and using 
Eqs. (58) and (67), we obtain the energy equation 

2g(mo5'o + E^iVo) + Si - a{n + 5) 



or equivalently the exphcit expression of the energy eigenvalues 



i;.(D) TlnlVo / ml /^ni V 1« 2 2 ^ 2 , 2\ 

_ 4(Vo^ - Sj) + /s:^, - SqmoSo 
"^"^ 4T/o^ + 

Kni = ga(2n + 1) ± ^Jq^a^ + 4{Sl - V^) + qa^{D + 2/ - 1)(D + 2/ - 3), 

(n = 0,l,2,3,---, Z = 0,l,2,---), (69) 

and 

q^a^ + qa^{D + 21- 2f + 45"^ > qa^ + AV^, (70) 

is a constraint over the potential parameters. In the case of pure vector potential {Sq — 

0, Vo 7^ 0), the energy equation (69) reduces to 



Kni = qa{2n + 1) ± q^a^ + qa^{D + 2/ - 1)(D + 2/ - 3) - W^, D>1, 

(n = 0,l,2,3,---, Z = 0,l,2,-..), (71) 

with the constraint 

q^a'^ + qa\D + 2l-2f>qa^ + 4V^, (72) 

for real bound states. Prom the above result, it is not difficult to conclude that the two energy 
solutions are valid for the particle and the second one corresponds to the anti-particle energy. 
The restriction that gives the critical coupling value leads to the result 



n<^\^^Aq^ml-V,'-^^-^[q+iD + 2l-l)iD + 2l-3)]-Vij-^, (73) 

i.e., there are only finitly many eigenvalues. In order that at least one level might exist, it 

is necessary that the inequality: 

qa + ^Jqa'' [q + {D + 2l- 1){D + 21- 3)] - AV^ < 2^Aq^ml-V^, (74) 
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is fulfilled. As can be seen from Eq. (73), there are at most only two lower- lying states 
{n — 0, 1) for the Klein-Gordon particle of mass unity when the parameter a — 1 and q — ±1 
for any arbitrary value of Vq. Therefore, we have the inequality: 



n < ± 



D + 21-2 



1 

2' 



(75) 



For instance, if one selects Vq = (-D + 2/ — 2)/2, then the real bound states must be restricted 
by n < (^^yl6-{D + 2l-2y - /2. 

Having solved the D-dimensional Klein-Gordon equation for scalar and vector general 
Hulthen-type potentials, we should make some remarks. 

(i) For s-wave {I — 0), the exact energy eigenvalues of the ID Klein-Gordon equation 
becomes 



2q 4q{W,' + Kl) 



[kI + 4 (V,^ - Si)] [{2So + Mmof - kI- Wi\ , (76) 



where 



= qa{2n + 1) + ^ q^a'' ^ 4{Sl -V^), 
_ 4(n' - Si) + kI- SqmoSo 

In order that at least one level might exist, it is necessary that the inequality 

16q'ml > r,l {AV^ + kI) , q^a" + ASl > AV^ 
is fulfilled. In the case of pure vector potential {Sq — 0,Vo 0), the energy spectrum 



(77) 



(78) 



" 2q 



with 



qa 



{2n+l) + y/q^a^-AVi 



rrin 



1 



^ AV^ + \qa{2n + 1) + ^/q''a^ - AV^ 



Wq'^ml > AVq + 



qa{2n + 1) + \J q^a^ - AV^ 



-n 2 



, qoi > 2Vo. 



16g2' 
(79) 

(80) 



We notice that the result given in Eq. (79) is identical to Eq. (31) of Ref. [30]. There are 
only two lower-lying states (n = 0, 1) for the Klein-Gordon particle of a rest mass mo = 1.0 
and a — 1.0. For example, one may calculate the ground state energy for the coupling 
strength Vq — qa/2 as 



En 



2q 



1 ± 



2mog 



(81) 



14 



Further, in the case of pure scalar potential {Vq — 0, Sq 0), the energy spectrum 

1 



4q 



qa {2n + 1) + y/q^a^ + 43^ 



qa{2n + 1) + ^Jq^a^ + ASi 



X i {2So + 4:moqf 



qa{2n + 1) + q^a^ + AS^ 



(82) 



In the above equation, all bound states appear in pairs, with energies ±i?„. Since the Klein- 
Gordon equation is independent of the sign of En for scalar potentials, the wavefunctions 
become the same for both energy values. We notice that Eq. (82) is identical to Eq. (24) of 
Ref. [47] obtained by N-U method and to Eq. (20) of Ref. [48] obtained by supersymmetric 
method. If the range parameter a is chosen to be o; = 1/Ac, where Ac = h/rriQC — 1/mo 
denotes the Compton wavelength of the Klein-Gordon particle. It can be seen easily that 
while ^0 — >^ in ground state (n = 0), all energy eigenvalues tend to the value Eq pa 0.866 
mo. 

(ii) For D — the mixed scalar and vector Hulthen potentials (g = 1, Z 7^ 0), the energy 
eigenvalues are 



E.. 



nl 



where 



-± 



a 



{2n + 1) + Vq;2 + AB 



AV^ + [a{2n + 1) + ^Ja^ + AB\ ' 



- ) , (83) 



4(K' - -^o) + [«(2^ + 1) + Va;2 + AB\ - SmoSo 



, B^ S^-V;^ + aH{l + l). (84) 



AV^ + [a{2n + 1) + Va^ + AB\ ' 
In order that at least one level might exist, it is necessary that the inequality 

16m^ > r]li i^AV^ + [a(2n + 1) + Va;2 + 4s]'| , a" + 4S > 0, 

is fulfilled. In the case of pure vector (5*0 — 0,Vo 0),the energy eigenvalues are 

Vo 



(85) 



a{2n + 1) + ^J a^{2l + 1)^ - AVi 



X 



rrin 



^ 4Vo2 + a{2n + 1) + ^/^^i2^TW^^^ 
with the following constraint over the potential parameters: 



16' 



(86) 



16m^ > AV^ + 



a{2n + 1) + ^Ja'^{2l + l)^ -AV^ 
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{21 + l)a > 2Vo. 



(87) 



A preliminary analysis about the possibility of obtaining real bound states for the Klein- 
Gordon particle of mass unity is done for various dimensions as follows: when D = 3, we 
have only two lower- lying real bound states {n — 0, 1) for the angular quantum number 
I — 0, one lower-lying state {n = 0) for Z = 1 and no states for Z > 2. When D = 4, we have 
two lower-lying states for Z = and no states for I > 1. Further, when D = 5, we have one 
lower-lying state for i = and no states for i > 1. Finally, when D > 6, there are no real 
bound states for / > 0, however, only scattering states will be possible. 

(iii) When D = 3 and / = 0, the centrifugal term (-P+^^~iK^^+2^~3) _ g^j^j ^j-^g g^p_ 
proximation term (^+^'~i^K^+^'^3)a e — _ Thus, letting / — and D = 3 in Eq. 

(10), it reduces to the exact spectrum formula and normalized radial eigenfunctions of the 
Klein-Gordon equation for vector and scalar general Hulthen-type potentials: 

/ 2 F2 _ 2gro(mo5o + EniVp) + ro{S^ - Vp^) _ n + 6 
V'"o 2q^{n + S) 9.r. ' 



2ro 



, (n = 0,l,2,3,---) 



(88) 



which gives 



K = ^ + (^qa{2n+l) + ^q^a^ + 4(5^ - Vi)j 



rrin 



\^ 4Vi + (qa{2n + 1) + V?'"' + 4(-S'2 - l^o')) 



16q^ml > r]^ 



Vn 



AV^ + (qa{2n + 1) + q^a'' + A{Sl -V^ 



|) , (n = 0, 1,2,3,...) 



gV + 4(5o'-Vo')>0, 



4(1^' - -S-o') + (g«(2n + 1) + VgV + 4(^2_y^2)^' _ g^^^^^ 



(89) 



+ (qa{2n + 1) + ^q^a' + 4(^2 - Vi)^ 

(iv) In the case — Vq Hulthen potential, Eqs. (67) and (68) can be reduced to the 
relativistic energy equation: 



^ ^ roH(^o + ^i?) ri + 5 . (D + 2/ - 1) 

"^0 - -^R = , fN tt;—^ (> = n 



(n + S) 

which is Eq. (22) of Ref. [19]. 



2ro 



, (n,Z = 0,l,2,...). (90) 
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(v) We discuss non-relativistic limit of the energy equation (90). When Vq — Sq, Eq. 
(10) reduces to a Schrodinger-hke equation for the potential 2V{r). In other words, the 
non-relativistic limit is the Schrodinger equation for the potential —2Voe~'^^'^° / [l — e~^/'^°] . 
Hence, using the transformation niQ + Er — 2mo and mo — Er — > —Ej^ji [22], we obtain 
the non-relativistic energy equation: 



Enr — 



57710 o; 



{n + 6) 



a 



(91) 



which is Eq. (23) of Ref. [19] with S is given in Eq. (90). It is worthwhile to remark that 
Eq. (91) is identical to Eq. (59) of Ref. [30] when the potential is 2V{r) and a becomes 
pure imaginary, i.e., a — > ia. 

Thus, in the weak couphng condition, [(n -|- 6)/morof <^ 1, [Voro/ {n + 6)]^ -C 1, expand- 
ing the energy equation (90), retaining only the term containing the power of (1/moro)^ and 
(?^oVb)^, we have the relativistic energy 

VoTo 



Er Enr + mo + 4mo 



(92) 



which is Eq. (24) of Ref. [19] with 5 is given in Eq. (90). The first term is the non-relativistic 
energy and third term is the relativistic approximation to energy. 

(vi) For a more specific case where q — —l, the usual Hulthen potential is reduced to the 
shifted usual Woods-Saxon (WS) potential 



V{r) = -Vo + 



Vo 



-, S{r)^-So + 



So 



-, ro = a 



1 + e-^'Ao ' ^ ' 1 + e-''/''o '' 

and hence the energy eigenvalues, Eq. (69), for the general WS-type potentials, i.e., q 
and then q = 1, become 

CnlVo 



(93) 

-q 



E. 



nl 



± 



^2a^ + 4:{Si - V^) - a'^{D + 21- 2)2 - a{2n + 1) 



rrin 



^ 4V^ + V2a2 + 4(^2 _ y2^ _ c^2p + 2/ - 2)2 - a{2n + 1) ^ 

IQml > Ci |4V^o' + ^2a2 + 4(^2 - Vq^) -a^{D + 2l- 2f -a{2n+ 1) 
where 

4(Vo' - -^2) + [V2«' + ^{Sl - Vi) - a\D + 21- 2f - a{2n + I] 



inl 

4 



(94) 



8mo5'o 



Wi+ y2Q;2 + ^{Sl - V^) - a^{D + 21- 2)2 - a{2n + 1) 



(95) 
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and the inequality 

2o? + A^Sl > AV^ + a^{D + 21- 2f, (96) 
must be fulfilled. In the pure vector potential [Sq = 0, Vq 7^ 0), the energy eigenvalues are 



E, 



(D) 
nl 



2 



^J2a^ - a^{D + 21- 2)2 - - a{2n + 1) 



X 



rrin 



^ + [V2«^ -a^{D + 2l- 2)2 - AV^ - a{2n + 1) 



16' 



16m5 > AV^ + 



^2a2 -a^^D + 2l- 2)2 - AV^ - a{2n + 1) 



where 



(97) 
(98) 



2a'^ >AV^ + a'^{D + 2l -2f. 

^From the above equation, for any given a, the Klein-Gordon equation with the usual shifted 
WS potential has negartive eigenvalues, i.e, E^i < 0. The given restriction in (99) imposes 
the critical coupling value and thus leads to the result 



n<- 

a 



(99) 



i.e., there are only finitly many eigenvalues. Further, it is necessary that the inequality: 



^a2 [2-{D + 2l-2Y]-AVi >a + 2^ Ami - K', 



(100) 



must be fulfilled. 

Now, let us find the wave function yn{s), which is the polynomial solution of 
hypergeometric-type equation, we multiply Eq. (19) by the weight function p{s) so that 
it can be rewritten in self-adjoint form [36] 



Hs)y'^{s)]' + Xp{s)yn{s)^0. 
The weight function p(s) which satisfies Eqs. (24) and (102) has the form 

p(s) =s2^"'(l-s)Vl+4{/?2+7)^ 

and consequently from the Rodrigues relation (23), we obtain 



(101) 



(102) 



Vnl 
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= S.,pf""^^^^^^\l-2.). (103) 

On the other hand, inserting the values of (7(s),7r(s) and t(s) given in Eqs. (59), (63) and 
(64) into Eq. (20), one can find the other part of the wave function as 

0(S) = S^-(l-s)M^+VW+7)]_ ^;L04) 

Hence, the wave function in Eq. (18) becomes 

9{s) = C„,5^-(l - .)^hV^^^(^]pf -'^^+^^^1 - 2s) 

= Cr,i/^\l - s)^Pi'^-''''-'V - 2s), s e [g,0). (105) 
Finally, the radial wave functions of the Klein-Gordon equation are obtained as 



where Nni is the radial normalization factor. Thus, the radial wave functions for the s-wave 
Klein-Gordon equation with pure vector Hulthen potential in ID reduces to 

Rn{r) = Cnie-^f<^^' (1 - ^e-^/^o)(<'+«)/(2?)pr°^"^°-''"""/'\l - 2qe-^l^^), (107) 

where a — yjq^ — ^Vq / which is identical to Eq. (35) of Ref. [30]. Finally, from Eqs. (5) 
and (6), the total wave functions for the usual Hulthen potential is 

= Af„ir-<°-"''e-V^<^(l - e-'-/'»)«pi"°^"'°-'=i'-''-"(l - 2e--/'«) 



/ (2n,,-i + 2m^+l)n,,-i! p{m' ,m'),^^a ^ nnsl 

V 2r (nz.-! + 2mO ^'^^^^"^^ (cos^^_i), (108) 

where 

<^ = ^ (l + + 2/ - 2)2 + ^rl{Sl - Vi)\^ , (/ = 0,1,2,---). (109) 

V. CONCLUSIONS 

We have analytically found an approximate bound state energy eigenvalues and their 
corresponding wave functions of the £)-dimensional Klein-Gordon equation for the spin- 
zero particle with the scalar and vector general Hulthen-type potentials using the N-U 
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method. The analytic energy equation and the wave functions expressed in terms of Jacobi 
polynomials and can be reduced to their well-known ID and 3D solutions. The relativistic 
energies E^i in £)-dimensions defined explicitly in Eq. (68) is for the particle and anti- 
particle energies for a given dimension D, quantum numbers n and / and also coupling 
constants satisfying the given particular constraints. For pure attractive scalar potential, 
all bound states appear in pairs, with energies ±Eni- Since the Klein-Gordon equation is 
independent of the sign of Eni for scalar potential, the wave functions become the same 
for both energy values. When I — 0, the results in this work reduce to exact solution of 
bound states of s-wave Klein-Gordon equation with vector and scalar general Hulthen-type 
potentials. A preliminary analysis about the possibility of obtaining finite number states for 
the Klein-Gordon particle of mass unity shows that the bound state energies are positive for 
the general Hulthen-type potentials and negative for general Woods-Saxon-type potentials. 

In Table 1, we have obtained numerical results for the binding energies of ground state 
for mass unity ID Klein Gordon equation with pure vector and pure scalar cases. In pure 
vector potential, for fixed coupling constants Vq = 0.25 mo and various values oi a — 0.5, 1.0 
and 2.0, most of the binding energies Eq are decreasing with increasing deformation constant 
q. Also the values obtained with pure scalar potential, 5*0 = 0.25 mo, are also given in Table 
1 for comparison. Furthermore, in Table 2, we give the binding energies for the n — 1, 
a — 0.5, 1.0 and n — 2,a — 0.5. Obviously, we have no bound states for n = 1 when a — 2.0 
and for n = 2 when a — 1.0 and 2.0. 
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TABLE I: Relativistic ground state binding energies, Eq, for a particle of mass, niQ = 1 given by 

Eqs. (79) and (82) for pure vector and pure scalar potentials, respectively, as a function of q for 

various values of a. 

Pure vector (Fq = 0.25) Pure Scalar {So = 0.25) 



Q 


a = 0.0 


a = 1.0 


^ on 
a = z.U 


a = U.O 


^ in 
a = l.U 


^ on 
a = z.U 


0.1 








0.755639 


0.947484 


0.946410 


0.5 




0.911438 


0.500000 


0.929812 


0.996314 


0.645297 


1.0 


0.820971 


0.970804 


0.250000 


0.986425 


0.964541 


0.480683 


1.5 


0.991673 


0.940945 


0.166683 


0.998606 


0.941246 


0.398341 


2.0 


0.999840 


0.923893 


0.125000 


0.999903 


0.926256 


0.347332 


2.5 


0.999140 


0.913089 


0.100000 


0.998350 


0.916087 


0.311871 


5.0 


0.988667 


0.890301 


0.050000 


0.988537 


0.892966 


0.222136 


7.5 


0.982893 


0.882371 


0.033334 


0.982996 


0.884414 


0.181787 


10 


0.979613 


0.878345 


0.025000 


0.979771 


0.879977 


0.157607 



'For real bound states q^a^>0.25. 
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TABLE II: Relativistic binding energies of the excited states, Ei and E2, for a particle of mass, 
mo = 1 given by Eqs. (79) and (82) for pure vector and pure scalar potentials, respectively, as a 
function of q for various values of a. 



n = 1 


n = 2 


n 


1 b 
= 1 




n = 2 


q a = 0.5 


a = l.l) a = u.o 


a 


= 0.0 


« in 
a = l.U 


a = U.O 


U.l — 




0. 


,995674 


0.771938 


0.922413 


U.o — 


U.83U948 - 


U, 


,984202 


0.5718/3 


0. 829156 


1.0 0.996421 


0.347292 0.880588 


0. 


,954903 


0.450227 


0.779514 


1.5 0.949420 


0.229259 0.769589 


0. 


,935491 


0.381304 


0.752337 


2.0 0.928534 


0.171407 0.737131 


0. 


,922604 


0.336188 


0.735135 


2.5 0.916027 


0.136934 0.719697 


0. 


,913600 


0.303882 


0.723317 


5.0 0.891025 


0.068342 0.688449 


0. 


,892282 


0.219316 


0.695608 


7.5 0.882692 


0.045546 0.678994 


0. 


,884103 


0.180255 


0.684996 


10 0.878525 


0.034156 0.674438 


0. 


,879801 


0.156613 


0.679406 



"Pure vector case with coupling constant Vb = 0.25. 
^Pure scalar case with coupling constant = 0.25. 
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